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^5 (54) Title: METHOD AND ARRANGEMENT FOR DETECTING AND MEASURING THE PHASE OF PERIODICAL BIOSIG- 
§H NALS 

^= (54) Bczeichnung: VERFAHREN UND ANORDNUNG ZUR DETEKTION UND ZUR MESSUNG DER PHASE VON PERIO- 
^ DISCHEN BIOSIGNALEN 

f 57 ) Abstract: The invention relates to a method and an arrangement for detecting and measuring the phase of periodical biosignals. 

O The aim of the invention is to detect and measure a causal phase response in periodical biosignals, with improved reliability and 
^ higher speed compared to conventional methods, simultaneously reducing the computing power required. According to the invention, 
a state observer is set up in parallel with the analysed biological system. The output variables of the biological system and the 
^ observer are evaluated in order to minimise the occurring error. The determined phase of a periodical biosignal can, for example, be 
O used for functional diagnostic purposes. 

§ (57) Zusammenfassung: Die Erflndung betrifft ein Verfahren und eine Anordnung zur Detektion und zurMessong der Phase von 
periodischen Biosignalen. Dabei ist es Ziel, einen kausalen Phasengang in periodischen Biosignalen mil einer gegenuber herkomm- 
lichen Verfahren besseren Zuveriassigkeit und hdheren Geschwindigkeitbei gleichzeilig reduziertem Rechenaufwand zu detekueren 
Q und zu messen. Erfindungsgemass wird parallel zum analysierten biologischen System ein Zustandsbeobachter aofgestellL Die 
^ Ausgangsgrossen des biologischen Systems und des Beobachters werden ausgewertet, am den auftrelenden Fehler zu minimieren. 
Die ermittelte Phase eines periodischen Biosignals kann beispielsweise ftir funktionsdiagnostische Zwecke genulzt werden. 
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Verfahren gnd Anordmang znr Detektion und zur Messung 
dar Phase von pexiodischen Biosignalen 



5 Die Erfindung betrifft ein Verfahren und eine Anordnung zur 
echtzeitf ahigen zuverlassigen Detektion und Messung der Phase 
periodischer physiologischer Grdfien Oder Biosignale. 

Im Stand der Technik sind Verfahren bekannt, die zur Bestim- 
10 mung der Phase ein uber den zeitlichen Verlauf des Biosignals 
gleitendes oder sequentiell angeordnetes Analysef enster 
verwenden. Auf den im Analysef enster liegenden Signalab- 
schnitt werden auf der Fourier-Transformation basierende 

Methoden oder deskriptive Statistiken angewandt. So werden 

15 beispielsweise in der Perimetrie periodisch auf leuchtende 
Lichtmarken definierter Intensitat und hinreichend hoher 
Frequenz (Qber etwa 4Hz) benutzt, um die Funktionsf Mhigkeit 
des visuellen Systems zu QberprOfen. Fur den Funktionstest 
wird das Elektroenzephalogramm (EEG) erfasst und die Reizant- 
20 wort auf den visuellen Stimulus bezOglich der Amplituden und 
der Phase analysiert. Die Phase der Reizantwort ist einer der 
entscheidenden diagnostischen Parameter in der Funktions- 
diagnostik. 

25 Bei den bisherigen Verfahren ist nachteilig, dass die statis- 
tische Unsicherheit der Detektion bzw. die Ungenauigkeit der 
Messung sehr hoch ist. Die Unsicherheit und die Ungenauigkeit 
ergeben sich aus der Signaltheorie als Folge von und im 
Zusammenhang mit der Lange des Analysef ensters . Die Theorie 

30 besagt f dass mit abnehmender LSnge des Analysef ensters die 
statistische Unsicherheit und damit die Ungenauigkeit zuneh- 
men, was in der praktischen Signalanalyse auch hinreichend 
bewiesen und bekannt ist. Fur ein statistisch besseres Ergeb- 
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nis miisste zunachst die FensterlSnge vergrfcftert werden. Aus 
der Physiologie jedoch ist bekannt, dass sich die Phase rela- 
tiv schnell andern kann und diese Anderungen auch diagnos- 
tisch relevant sind. Bei einem langen Analysefenster geht die 
5 wertvolle Information uber die Phasenanderung verloren und 
die statistische Unsicherheit des Messergebnisses nimmt 
infolge der Anderungen nicht zwangslaufig ab. 

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren und 
10 eine Anordnung anzugeben, mit denen es moglich ist, einen 
kausalen Phasengang in periodischen Biosignalen mit einer 
gegenuber herkGmmlichen Verfahren besseren Zuveriassigkeit 
und h&heren Geschwindigkeit bei gleichzeitig reduziertem 
Rechenaufwand zu detektieren und zu messen. 

15 

Erf indungsgemaft gelingt die Losung der Aufgabe dadurch, dass 
periodische Biosignale entsprechend ihrem physikalischen und 
physiologischen Ursprung erfasst werden, dass parallel zum 
analysierten biologischen System ein Zustandsbeobachter 
20 aufgestellt wird, dass die AusgangsgrOfcen des biologischen 
Systems und des Beobachters mit Hilfe eines Kalman-Filters 
ausgewertet werden und zur Bestimmung der Phase genutzt 
werden. 

25 Bei dem erf indungsgemaJJen Verfahren wird die Phase eines 
periodischen Biosignals ermittelt und fUr funktionsdiagnosti- 
sche Zwecke genutzt. So ist beispielsweise eine verlangerte 
Phase im Vergleich mit der gesunder Versuchsobjekte ein wich- 
tiger Hinweis auf funktionelle Probleme des untersuchten 

30 biologischen Systems. 

Bei der erf indungsgemaften Anordnung ist parallel zum unter- 
suchten biologischen System, das mit einem Zustandsmodell 
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nachgebildet wird, ein Zustandsbeobachter angeordnet, der . 
entsprechend dem Systemmodell die Zustandsgroiie Phase basie- 
rend auf dem Kalman-Filter schatzt. 

Von Vorteil ist hierbei, dass die Schatzung der Phase 
kontinuierlich stattfinden kann und kein gleitendes Oder 
sequentiell angewandtes Analysef enster notwendig ist. Damit 
wird die Analyse , der zeitlichen Phasenanderungen erst 
m6glich. Im Gegensatz zum relativ komplizierten theoretischen 
Hintergrund dieses Phasenschatzers ist die praktische Umset- 
zung einfach. Im Vergleich zu herkSmmlichen Verfahren bend- 
tigt sie wesentlich weniger Rechenleistung, so dass eine 
Phasenschatzung in Echtzeit mttglich ist. 

Die Erfindung wird im folgenden anhand der theoret-i-SGhen- 
Herleitung und eines Ausf iihrungsbeispieles naher erlautert. 
In den zugehorigen Zeichnungen zeigen: 

ein Blockdiagramm eines Beobachterkonzeptes; 

ein Systemmodell einer erf indungsgemafien Anordnung; 

ein prinzipielle Darstellung des Zustandsbeobach- 
ters zur Messung der Phase in periodischen Biosig- 
nalen; 

einen Verlauf einer geschatzten Phase far eine 
Harmonische der Frequenz 8Hz und. Phase 2rad fUr die 
statischen Kalman-Faktoren 2 und 20; 

einen Verlauf einer geschatzten Phase fUr eine 
verrauschte Harmonische der Frequenz 8HZ und Phase 



Fig. 1 
Fig. 2 
Fig . 3 

Fig. 4 

Fig . t 5 
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2rad mit einem SNR von OdB (unten) und dynamischem 
Kalman-Faktor (oben) ; 

Fig. 6 eine Phasenschatzung des Signals wie in Fig. 5 mit 
5 statischen Kalman-Faktoren; 

Fig. 7 Ergebnisse der Phasenschatzung (rechts) an realen 
Signalen (links) ; 

10 Fig- 8 eine additive Oberlagerung einer Harmonischen einer 
Frequenz von 8Hz mit dem Rauschen und gezielter 
Verstimmung der Analysef requenz (oben) und den 
Phasenverlauf mit Anstieg (unten) . 

15 

Ein biologische System, das ein periodisches Biosignal produ- 
ziert bzw. auf ein periodisches Eingangssignal antwortet, ist 
in Fig. 1 als Zustandsmodell „reales System" dargestellt. Die 
nachfolgend genannten Zustandsgleichungen (1) und (2) 
20 beschreiben dieses System (Fettgedruckte grolie Buchstaben 
stehen f iir Matrizen, kleine far Vektoren) : 

x(t) = A-x(t) + Bu(t); x(0) = x(t 0 ) (1) 
y(t) = C-x(t). (2) 

25 

Fur die weiteren Betrachtungen wird ein additives Signalmo- 
dell angenomraen, welches eine harmonische Schwingung und 
weiftes normal verteiltes Rauschen summiert: 

30 y(t) = ysin((ot + (p(t)) + r p (t) (3) 

Das Ziel besteht darin, ein Systemmodell zu konstruieren, 
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dessen Variable x(t) die Phase <p(t) des zu untersuchenden 
Signals y(t) reprasentiert . Die Phase kann nicht direkt 
gemessen werden, da sie Argument einer trigonometrischen 
Funktion ist. Daher wird eine Hilf skonstruktion benotigt. 
5 Eine solche Konstruktion ist ein Zustandsbeobachter , der 
parallel zum untersuchten System angeordnet ist. Der Beobach- 
ter schatzt die Zustandsvariable mit Hilfe der Minimierung 
einer Fehlerfunktion, welche die Ausgange des realen Systems 
und des Beobachters vergleicht. Auf diese Weise kann nach 
10 abgeschlossener Fehlerminimierung die Zustandsvariable Phase 
direkt gemessen werden. 

Das Blockdiagramm des Beobachterkonzeptes ist in Fig. 1 
dargestellt. Da x(t) nicht direkt gemessen werden kann, wird 
15 3Cm(t) im Beobachter geschatzt. Die innere Schleife im 
Beobachter minimiert. den Fehler von y n (t) beziiglich y(t) mit 
Hilfe der Korrekturmatrix K. FUr den Beobachter ergeben sich 
dann die Zustandsgleichungen (4) und (5): 

20 x M (t) = A • x M (t) + B-u(t) + K- [y(t) - y M (t)) , (4 ) 

y M (t) = Cx M (t). (5) 

Aus (4) und (5) folgt: 

25 x M (t) = (A-KC).x M (t) + B.u(t) + K-y(t). (6) 

Es wird davon ausgegangen, dass beide Systeme unterschiedli- 
che Anf angsbedingungen haben. Daraus ergibt sich der Beobach- 
tungsfehler: 

30 

e(t) = x(t)-x M (t). (7) 
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Der Beobachtungsfehler verschwindet iterativ mit Hilfe der 
Korrekturmatrix K, so dass 

e (t) = 0 for t->oo. (8) 

Die Dynamik und die Stabilitat der Schatzung konnen mit der 
Dif ferentialgleichung des Beobachtungsf ehlers (9) beschrieben 
werden: 



e(t) = x(t)-x M (t). O) 

Durch Umstellung und weitere Zwischenschritte erhalt man: 

15 c(t) = (A-KC)e(t). (10) 

Entsprechend dem Signalmodell (3) ist damit zu rechnen, dass 
das untersuchte Signal durch Rauschen gestdrt ist. Urn den 
Einfluss des Rauschens zu reduzieren, wird ein Kalman- Filter 
20 eingesetzt. Das, Rauschen berttcksichtigend wird das System 
durch folgende Zustandsgleichungen beschrieben: 

Syst emzust and : x(t) = A • x(t) + B • u(t) + M • r s (t) ( 1 1 ) 

Systemausgang: y(t) = Cx(t) + r p (t) (12) 
2 5 Beobacht er : x M (t) = A • x M (t) + B - u(t) + K(t) ■ [y(t) - C • x M (t)] , (13) 

wobei 

K(t) ist die korrekturmatrix, welche zu erreichen hat, 
30 dass e(t) = x(t)-x M (t)-*0, 

o(t) ist der Beobachtungsfehler, 
x s (t)ist das Systemrauschen, und 
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r p (t) ist das Prozessrauschen. 



Um die Herleitung zu vereinf achen, wird angenommen, dass die 
5 Rauschanteile breitbandige GaufJsche Nullmittelprozesse mit 
bekannten Kovarianzen sind: 

cov fp (t ls t 2 ) = Efcft) • r p T (t 2 )}= R p (Q - 8(t, - 1 2 ) 
cov Jt,, t 2 ) = E{r 8 (t,) -r s T (t 2 )} = R,(t,) .5(t, - 1 2 ) 

10 Die Rauschanteile sind voneinander unabhangig, also 

cov fp( .(t M t 2 ) = 0. (15) 

Far eine konsistente Schatzung von x(t) muss die Fehlerleis- 
15 tung mit Hilfe der Matrix K(t) minimiert werden: 

E{e T (t) ■ e(t)} = E{ef (t) + e*(t) + • + eJS (t)}= f (K(t))iMin. , (16) 

Unter Annahme der stochastischen Beziehungen bezuglich der 
20 Kovarianzen wird entsprechend dem Kalman-Filter eine geeig- 
nete Korrekturraatrix K(t) abgeleitet: 

K(t) = cov e (t).C T .R; ! (t). (17) 

25 Die Formel ftir die Fehlerkovarianz cov E (t) kann vom Kalman- 
Filter abgeleitet werden: 

covc(t) = Acov c (t) + cov c (t)^ (18) 
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Schatzung der Phase: 

Das untersuchte Signal wird entsprechend (3) aus der Summe 
einer Harmonischen und Rauschen modelliert: 

5 

y(t) = yjt) + r p (t) = y - sin(cot + <p(t))+ r p (t),co = 2nf . (19) 

Die Phase ergibt sich aus der Dif f erentialgleichung : 

10 9(0 = -a(p(t) + r s (t) a>0. (20) 

Daraus ergibt sich das Systemmodel in Figur 2. Die Phase kann 
nicht direkt gemessen werden. Daher liegt parallel zum System 
ein Beobachter, in dem ein direkter Zugriff auf die 

15 geschatzte Phase cp M (t) moglich ist. Allerdings der nichtline- 
are Anteil y n i(t) in (19) ist ftir das Bepbachterkonzept 
ungiinstig. Eine geeignete Linearisierung yi(9(t) r t) wird 
benotigt, um eine linear e Beziehung zwischen der Zustandsva- 
riablen q> M (t) und dem Ausgang y M (t) entsprechend (5) herzu- 

20 stellen. Basierend auf der Taylor-Linearisierung kann der 
Beobachter f olgendermafien f ormuliert werden: 

<i> M (0 = -a ■ 9m (0 + K(t) - (y (t) - y M (t)) (21) 

25 y M (t) = yi(<P M (0,t). (22) 

Entsprechend (21) und (22) wird der Beobachter modelliert:, 
wie in Fig. 2 dargestellt . Im Ergebnis der Linearisierung im 
Arbeit spunkt q> B , wird (22) mit (5) verkniipft. Daraus ergibt 
30 sich der Faktor C, welcher in (17) verwendet wird, um den 
Korrekturfaktor K(t) zu bestimmen: 



BESTATIGUNGSKOPIE 



\yO 2004/04761 




PCT/EP2003/013355 



- 9 - 



K(t) = cov e (t) ■ y ■ cos(o)t + q> 9 (t)) • R; ! (t). 



(23) 



Als Arbeitspunkt wird die zu be?stimmende Phase gewahlt 



5 



q> B ( t ) = q> M ( t )> 



(24) 



und die resultierende Dif f erentialgleichung fur die Phase ist 
ci>M (0 = -a ■ <P M (0 + cov c (t) • y • cos(cot + <p M (t)) - [y(t) - y • sin(cot + <p M (t))]- Rp 1 (t). ( 25 ) 

10 

Entsprechend (25) kann der Phasenschatzer modelliert werden, 
wie in Fig. 3 dargestellt. 

Zur Phasensch&tzung in y(t) muss die Fehlerkovarianz berech- 
15 net werden. Aus (18) folgt: 

cove (t) = -2a • cov e (t) - y 2 • cos 2 (cot + <p M (t)) - cov 2 (t) - R; ! (t) + R 5 (t). (26) 

Gleichung (26) ergibt eine einfache Losung, falls hoherfre- 
20 quente Anteile in der Fehlerkovarianz nicht berQcksichtigt 
werden. Basierend auf (27) 



Bei geeigneter Wahl des Parameters a in (28) werden hochfre- 
30 quente Anteile in Folge temporaler Integration unterdrtickt, 
also es liegt ein Tief passverhalten vor. Unter Berticksichti- 



cos 2 (cot + cp M (t)) = 1 12 + cos(2(Dt + 2<p M (t)) 



(27) 



25 



kann (26) vereinfacht werden zu: 



covc(t) = -2a - cov c (t) - 1 12 • y 2 - R; f (t) -cov^(t) + R,(t) . 



(28) 
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gung des Tief passes kann (25) vereinfacht werden: 

<Pm (0 = -a (P M (t) + cov c (t) - y • cos(e)t + <p M (t)) • y(t) • R;'(t), (29) 

5 Damit wird der Beobachter, dargestellt in Fig. 3, verein- 
facht. Das in (29) vorgeschlagene System kann insbesondere 
zur Phasenschatzung von Harmonischen im Rauschen eingesetzt 
werden. 

10 In Fig. 4 ist der Verlauf der geschatzten Phase ftir eine 
Harmonische der Frequenz 8Hz und Phase 2rad f iir unterschied- 
liche Kalman-Faktoren dargestellt. Die Kalman-Faktoren sind 
statisch und betragen 2 bzw. 20. Wie der Grafik entnommen 
werden kann, mit einem niedrigeren Kalman-Faktor wird die 

15 Schatzung langsamer. Statische Kalman-Faktoren miissen dort 
eingesetzt werden, wo der Zeitpunkt der PhasenSnderung nicht 
bekannt ist. 

In Fig. 5 ist der Verlauf (unterer Teil) der geschatzten 
20 Phase fur eine verrauschte Harmonische der Frequenz 8Hz und 
Phase 2rad mit einem SNR (Signal-Rausch-Verhaltnis) von OdB 
und dynamischen Kalman-Faktor (oberer Teil) dargestellt. Ist 
der Zeitpunkt der PhasenSnderung bekannt, so kann der Kalman- 
Faktor so konstruiert werden, dass zunSchst der Anderung 
25 gefolgt und anschlieBend die Varianz der Schatzung minimiert 
wird. 

Zum Vergleich ist in Fig. 6 die Phasenschatzung des " selben 
Signals wie in Fig. 5 mit statischen Kalman-Faktoren darge- 
30 stellt. 
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In Fig. 7 sind die Ergebnisse einer Phasenschatzung (rechte 
Spalte der Grafik) an realen Signalen dargestellt. Zur Stimu- 
lation eines beispielhaft untersuchten visuellen Systems 
wurden Lichtpulsf olgen mit der Wiederholrate von 8 Pulsen pro 
5 Sekunde verwendet, wobei abwechselnd einer Ruhepause (Zeit 
von 0 bis 2 der Zeitverlaufe der Signale* oben links und unten 
links) die Lichtstimulation (Zeit 2 bis 5) folgte. Darge- 
stellt ist der Verlauf des EEG (Elektroenzephalogramm, linke 
Spalte der Grafik) von zwei okzipitalen Positionen nach 16- 
10 facher reizbezogener Mittelung. Beide Zeitverlaufe zeigen 
einen deutlichen Sprung der Phase nach Beginn der Lichtstimu- 
lation. 

Die Phasenschatzung wird problematisch bei stark verrauschten 
15 Signalen. Generell gilt zwar, dass die Phase robuster gegen 
Storungen ist als die Amplituden, wie dies schlieBlich auch 
in der Inf ormationstechnik bekannt ist. Allerdings ist in 
diesem Grenzbereich zunachst die Frage nach dem Vorhandensein 
- also der Detektion - einer kausalen Phase zu klSren, erst 
20 dann ware die Phase zu schStzen. 

In Fig. 8 ist eine Harmonische der Frequenz 8Hz additiv dem 
Rauschen tiberlagert beginnend bei t=4s, wobei das SNR -lOdB 
betragt (oberer Verlauf in der Grafik) . Im Zeitbereich ist 

25 die Harmonische beztiglich ihrer Amplituden nicht nachweisbar. 
Verwendet man den Phasenschatzer mit einer gezielten Verstim- 
mung, hier mit einer Frequenz von 7.8Hz, also iam 0.2Hz 
weniger als die Frequenz der Harmonischen, so ergibt sich im 
Falle einer kausalen Phase ein Anstieg von 1.2rad/s (unterer 

30 Verlauf der Grafik) . Dieser Anstieg kann in Kombination mit 
einem Diskriminator direkt zur Detektion des Signals genutzt 
werden. 
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Bezugs zeichenli s to 



5 

a - Systemparameter, im Beobachtermodell wahlbar 

A, B, C, K - Matrizen im Zustandsmodell eines Systems 

(p(t) - Phase im Systemmodell, zu schatzende Gr6fte 

q> M (t) - Phase im Beobachtermodell, messbare GroBe 

10 r p (t) - Prozessrauschen 

r s (t) - Systemrauschen 

u(t) - Eingangsvariable eines Systems im Zustandsmodell 

x(t) - Zustandsvariable eines Systems im Zustandsmodell 

xn(t) - Zustandsvariable des Beobachters im Zustandsmodell 

15 y(t) - Ausgangsvariable eines Systems im Zustandsmodell 

yk(t) ~ Ausgangsvariable des Beobachters .im Zustandsmodell 
yi(<p(t),t) - Linearisierungsoperator fttr Phase 
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Paten tanapriiclie 



5 
10 
15 
20 
25 

2. 

30 



Verfahren zur Detektion und zur Messung der Phase von 
Antwortsignalen (y(t)) eines Biosystems, folgende 
Schritte umfassend: 

a) Multiplikation des Antwortsignals (y(t)), dessen 
Phase (q>(t)) bestimmt werden soil, mit einem ersten 
Faktor; 

b) Multiplikation des aus Schritt a) erhaltenen Produkts 
mit einem zweiten Faktor, welcher durch eine trigono- 
metrische Funktion reprasentiert ist, deren Argument 
sich aus dem Produkt der Frequenz des untersuchten 
Antwortsignals mit der Zeit, addiert mit der gemesse- 
nen Phase, ergibt, wobei die Frequenz der trigonomet- 
rischen Analysef unktion der Frequenz entspricht, bei 
der die Phase ermittelt werden soli, oder urn einen 
bekannten Betrag von dieser Frequenz abweicht, 

c) Multiplikation der gemessene Phase mit einem dritten 
Faktor (a) ; 

d) Dif ferenzbildung aus dem im Schritt b) erhaltenen 
Produkt und dem im Schritt c) erhaltenen Produkt; 

e) Integration der im Schritt d) erhaltenen Dif ferenz 
Uber die Zeit, wobei das Ergebnis dieser Integration 
die zu ermittelnde Phase des Signals reprasentiert 

f) Wiederholung der Schritte a) bis e), bis ein Abbruch- 
kriterium erfiillt ist. 

Verfahren nach Anspruch 1, dadurch gekennzeichnet , dass 
der erste Faktor zeitlich konstant oder ver£nderlich 
gewShlt . wird. 4 
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3. Verfahren nach Anspruch 1 oder 2, dadurch qekennzeichnet , 
dass bei der Dif f erenzbildung im Schritt d) das Produkt 
aus Schritt c) von dem Produkt aus Schritt b) abgezogen 
wird. 

5 

4. Verfahren nach einem der Anspruche 1 bis 3, dadurch 
qekennzeichnet , dass das Antwortsignal (y(t) ) einem 
Zustandsbeobachter zugefuhrt wird, welcher die Verfah- 
rensschritte a) bis f) ausfuhrt, urn eine geschatzte Phase 

10 (q>n(t)) • zu bestimmen, wobei das Verfahren abgebrochen 

wird, wenn das Beobachterausgangsignal (y M (t)) weniger 
als einen durch eine Fehlerf unktion (cov e (t)) vorgegebe- 
nen Fehlerwert (e(t)) vom Antwortsignal (y(t) ) abweicht 
und wobei nach Abbruch des Verfahrens die geschatzte 

15 Phase (<p M (t)) gleich der Phase (<p(t)) des Antwortsignals 

gesetzt wird. 

5. Verfahren nach Anspruch A, dadurch gekennzeichnet , dass 
die Schritte a) bis d) im Beobachter nach der folgenden 

20 Formel ausgefiihrt werden: 

^,(0 = -a-P*(0+co v e(0^ . 

6. Anordnung zur Detektion und zur Messung der Phase von 
Antwortsignalen eines Biosystems r dadurch qekennzeichnet , 

25 dass die Anordnung einen Zustandsbeobachter umfasst, dem 

parallel zum Biosystem das Antwortsignal des untersuchten 
Systems eingespeist wird und welcher die Verfahrens- 
schritte nach einem der Anspruche 1 bis 5 ausfiihrt. * 

30 7. Anordnung. nach Anspruch 6, dadurch qekennzeichnet , dass 
der Zustandsbeobachter einen Kalman-Filter umfasst, mit 
welchem Storsignale aus dem Antwortsignal gefiltert 
werden . 
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